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ABSTRACT 

The paucity of observed dwarf galaxies in the Local Group relative to the abundance of predicted 
dark matter halos remains one of the greatest puzzles of the ACDM paradigm. Solving this puzzle 
now requires not only matching the numbers of objects but also understanding the details of their 
star formation histories. We present a summary of such histories derived from the HST data using 
the color-magnitude diagram fitting method. To reduce observational uncertainties, we condense the 
data into five cumulative parameters - the fractions of stellar mass formed in the last 1, 2, 5, and 10 
Gyr, and the mean stellar age. We interpret the new data with a phcnomcnological model based on 
the mass assembly histories of dark matter halos and the Schmidt law of star formation. The model 
correctly predicts the radial distribution of the dwarfs and the fractions of stars formed in the last 5 
and 10 Gyr. However, in order to be consistent with the observations, the model requires a significant 
amount of recent star formation in the last 2 Gyr. Within the framework of our model, this prolonged 
star formation can be achieved by adding a stochastic variation in the density threshold of the star 
formation law. The model results are not sensitive to late gas accretion, the slope of the Schmidt law, 
or the details of cosmic reionization. A few discrepancies still remain: our model typically predicts 
too large stellar masses, only a modest population of ultra-faint dwarfs, and a small number of dwarfs 
with anomalously young stellar populations. Nevertheless, the observed star formation histories of 
Local Group dwarfs are generally consistent the expected star formation in cold dark matter halos. 
Subject headings: cosmology: theory — galaxies: formation — galaxies: dwarf 



1. INTRODUCTION 

The discrepancy between the number of observed 
dwarf satellite galaxies of the Milky Way and the 
anticipated number of dark matter halos in cosmo- 
logical simulations ha s been heralded as the "miss- 
ing satellite problem" (iKlypin et ail 119991 : iMoore et al.1 
1999] iKauffmann et all 119931) and it still remains one 
of the greatest puzzles of the ACDM paradigm. One 
class of suggested solutions to this puzzle involves 
modifications to the nat ure of dark matter (e.g., 
ISpergel fc Steinhardtj l2000f ) or to the initial conditio ns 
of cosmic structure (e.g.. IZentner fc Bullockl l2003aH bl: 
iKamionkowski fe Liddld l2000t IColin et al.ll2000D . while 
another class of solutions invokes astrophysical argu- 
ments such as inefficient cooli ng of cosmic gas and/or 
feedback from youn g stars (e.g. lThoul fc W cinbcrj H996t 
iDekel fc Wol3l2003lh 

Motivated by this lat ter class of solutions, 

iKravtsov. Gnedin. fc Klvpinl (|2004l hereafter referred 
to as KGK04) developed a star formation model for 
the satellite galaxies, based on the mass assembly 
histories of individual dark matter halos in a ACDM 
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cosmological simulation. That model incorporated the 
accretion of gas in hierarchical mergers and the loss of 
gas caused by the extragalactic UV background, and 
included both a continuous mode and a starburst mode 
of star formation. The model correctly reproduced the 
observed number of satellites of the Milky Way and 
M31, the radial distribution of luminous satellites, and 
the morphological segregation of dwarf irregular (dlrr) 
and dwarf sph eroidal (dSph ) galax ies. Recent V-body 
simulations of iMadau et al.1 ()2008l ). with the updated 
WMAP3 cosmological parameters, have confirmed the 
qualitative picture of KGK04. 

The amount of observational data for Local Group 
dwarfs has been increasing steadily in the last several 
years, shifting the focus from simply counting the num- 
ber of dwarf galaxies to deriving their detailed star for- 
mation histories (SFHs). Most Local Group dwarfs 
have now been observed with the HST 7 which has pro- 
vided color-magnitude diagrams of resolved stellar pop- 
ulations in single or m ultiple fields (jDolphin et al.l [20051 : 
iHoltzman et al.ll2006l ). We parametrize these data as the 
fractions of stars formed in the last 1 Gyr, 2 Gyr, 5 Gyr, 
and 10 Gyr, and present these fractions in our Table [TJ 
The data show a great variety of star formation histo- 
ries - some continuous, some bursty, some truncated. In 
this paper we extend the model of KGK04 to account for 
such varied histories of the Local Group dwarfs and argue 
that adding a stochastic threshold to the star formation 
law greatly improves the comparison between models and 
observations. We also make predictions for the num- 
ber of ultra- faint satellites of the Milky Way and M31 
of which there have been many recent discoveries with 
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the Sloa n Digital Sky Survey ( S PSS) and the MegaCam 
survey (iBelokurov et all 12009. l2007t llrwin et aIT | 2007|: 
Liu et al.H2q08HMartin et al.ll2006t ISimon fc Gehall2007t 
Walsh et all 120071: iWillman et all l2005allhl : iZucker eTatl 
2004 l2006al fbl. 120071 k 

Our approach is complementary to running full hydro- 
dynamic simulations of galaxy formation. We attempt 
to model star formation with a phcnomcnological recipe, 
containing several but not too many free parameters. We 
investigate the effects of varying each of these parameters 
and constrain them by comparing the model predictions 
for each parameter set with the observations. In this way, 
we obtain empirical rules that govern star formation in 
dwarf galaxies, even though we may not yet fully under- 
stand their physical origin. In an alternative approach, 
one uses cosmological hydrodynamic simulations with a 
detailed prescription for gas cooling and star formation, 
and thus obtains more physically sound results. Dis- 
crepancies between the simulation results and the data 
may, however, be difficult to interpret or computation- 
ally expensive to investigate. By combining the two ap- 
proaches, our empirical rules will guide future develop- 
ments of more sophisticated models of galaxy formation. 

The paper is organized as follows. We give a brief sum- 
mary of the observational data in $2] We describe our 
model in Sj3] and report the primary results from adding 
stochasticity to the star formation law in SjT] In <JS]we de- 
scribe the results of several variants of our star formation 
model, including late accretion of gas, steeper slopes of 
the Schmidt law, an extended epoch of reionization, and 
rejection of galaxies with delayed star formation; none 
of these changes significantly alters our results. In Sj6]we 
show our model projections for the new low-mass dwarfs, 
and in JJ7]we summarize our main results. 

2. OBSERVED STAR FORMATION HISTORIES 

The star formation data used in this study, and listed 
in Table [TJ come primarily from Hubble Space Telescope 
observations with the WFPC2 camera. For a majority 
of the dwarfs, the SFHs were measured by t he color- 
magni tude diagram fitting method, described in|Dolphin 
(2002), using photometric d ata from the Local G roup 
Stellar Photometry Archive 1 (jHoltzman et al.ll2009) . We 
condense the data into five parameters - the fractions of 
stellar mass formed in the last 1, 2, 5, and 10 Gyr, re- 
spectively (/i G , f 2 c, fza, Aog), and the mass-weighted 
mean age of the stellar population (r) . The HST observa- 
tions included in archive are far from homogeneous, and, 
therefore, the uncertainties in the derived star formation 
histories cover a broad range. However, by concentrating 
on cumulative star formation fractions the uncertainties 
are greatly reduced, compared to specific star formation 
rates at specific times. Thus, the entries in Table [JJ do 
not have associated errors. Constraining comparisons 
between models and observations to cumulative distribu- 
tions greatly reduces the sensitivity of these comparisons 
to uncertainties in individual galaxies. 

Note that the observed HST fields do not necessarily 
cover the full extent of each galaxy and, therefore, do 
not allow us to calculate the total stellar mass. For the 
estimate of total stellar mass, with a few exce ptions, we 
take the values quoted in lDekel fc Wool (|2003l ). 

1 http:/ /astronomy.nmsu.edu/holtz/archival/html/lg.html 



The rest of the SFHs are taken from a variety of 
sources, as indicated in the notes to Table [TJ In the case 
of the LMC, the star form ation data are inferred from 
ISmecker-Hane et al.l (|2002l ). Since they give the star for- 
mation rates per unit area, in units of MQyr _1 deg _ , 
separately for the disk and the bar of the LMC, we as- 
sume that the extent of the disk is 3.5 times larger than 
the extent of the bar in determining the star formation 
fractions and the mean age (A. A. Cole, personal com- 
munication) . 

The star formation histories for five dwarfs (M32, M33, 
Antlia, Sextans, Urs a Minor) are taken from Table 1 of 
IDolphin et all (12001 . Since that table does not include 
the fraction of a galaxy's star formation in the last 2 Gyr 
or 5 Gyr, our table is also missing these entries, except for 
Sextans and Ursa Minor. IDolphin et al.l (|2005f ) reported 
that Sextans and Ursa Minor have formed no stars in 
the past 10 Gyr (i.e., /iog = 0), which implies that the 
fractions of stars formed in the last 1 Gyr, 2 Gyr, and 
5 Gyr are also consistent with zero. 

For the purpose of comparing the radial distributions 
of different morphological types with model predictions, 
we combine the Sc, Irr, and dlrr types into one broadly- 
defined "dlrr" group, and the dE, dSph types into the 
"dSph" group. 

One object not incl uded in Tabled] is th e Canis Major 
dwarf, discovered by iMartin et all (|2004f) in the SDSS 
field, which could be part of a larger Monoceros tidal 
stream. This object is the closest discovered to the 
Galaxy, at a heliocentric distance of only 8 kpc. Whether 
this galax y is still grayitatio nally self-bound is unclear at 
present ([Butler et al.l l2007h . but in any case it is being 
strongly tidally disrupted. An additional reason for not 
including it in our comparison is that such an object 
would have been completely disrupted in the ./V-body 
numerical simulation used to construct our models. 

Note also that Table [TJ does not include any of 
the newly-discovered SDSS and MegaCam ultrafaint 
dwarfs, since they are likely to have stellar masses be- 
low 5 x 10 5 M©. We list these new objects in Tabled] 
and compare their estimated stellar mass function with 
the predictions of our fiducial model in SjB] 

3. STAR FORMATION MODEL 

The KGK04 model of star formation in dwarf galax- 
ies is based on the mass assembly history of dark matter 
halos in a collisionless ACDM simulation of the Local 
Group- like environment. The simulation volume con- 
tains three large host halos, with virial masses (1.2 — 
1.7) x 10 12 hr 1 M Q at z = 0, resolved with - 10^ dark 
matter particles. The model incorporates the accretion 
of gas in hierarchical mergers and the loss of gas caused 
by the extragalactic UV backgroun d (following the filter- 
ing mass approach of[Gncdin 200j3), and includes both a 
continuous mode and a starburst mode of star formation. 

The model assumes that in the satellite halos the ac- 
creted gas dissipates its energy via radiative cooling and 
forms a disk. The surface density of the gas follows an 
exponential profile, 

= s o exp (-r/r d ), (1) 

with the scale length r d set by the satellite halo's virial 
radius r v - 1T and angular spin parameter A: 

r d = Ar vir 2- 1 / 2 exp [c^/K^) 2 ] . (2) 
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TABLE 1 

Star Formation Histories of Satellite Galaxies of MW and M31 



Galaxy 


Alternate Name 


Type 


Host 


fhost (kpc) log(M,, Q ) 


JIG 




J2G 


/5G 


/lOG 


T (Gyr) 


M33 


NGC598 


Sc 


M31 


203 


9.9 


0. 


093 










0. 


52 


8.4 


LMC 




lrr 


MW 


50 


9.7 


0. 


078 





.17 


0. 


42 


0. 


70 


6.7 


SMC 




Irr 


MW 


63 


9.2 


0. 


096 





.18 


0. 


.48 


0. 


65 


6.6 


M32 


NGC221 


dE 


M31 


6 


9.1 


0. 


042 










0. 


50 


8.5 


NGC205 


MHO 


dE 


M31 


58 


9.0 


0. 


0049 





.0050 


0. 


.055 


0. 


.48 


10.5 


IC10 




dlrr 


M31 


255 


8.7 


0. 


060 





.14 


0. 


52 


0. 


.75 


7.1 


NGC6822 




dlrr 


MW 


500 


8.7 


0. 


087 


0. 


.16 


0. 


57 


0. 


68 


6.9 


NGC3109 




dlrr 


MW 


1360 


8.7 


0. 


054 





.065 


0. 


084 


0. 


12 


11.1 


NGC185 




dSph 


M31 


175 


8.6 


0. 


0053 





.0053 


0. 


090 


0. 


51 


10.5 


NGC147 




dSph 


M31 


101 


8.4 


0. 


032 





.036 


0. 


050 


0. 


17 


12.4 


IC1613 




dlrr 


M31 


505 


8.3 


0. 


059 


0. 


.11 


0. 


42 


0. 


64 


7.7 


WLM 


DD0221 


dlrr 


M31 


840 


8.2 


0. 


14 





.35 


0. 


.55 


0. 


69 


6.7 


Sex B 


DDO70 


dlrr 


MW 


1320 


8.2 


0. 


049 





.067 


0. 


11 


0. 


21 


11.1 


Sex A 


DD075 


dlrr 


MW 


1440 


7.9 


0. 


15 





.29 


0. 


38 


0. 


11 


9.3 


Sagittarius 




dSph 


MW 


28 


7.7 


0. 


0008 





.0008 


0. 


52 


0. 


86 


6.5 


Fornax 




dSph 


MW 


138 


7.5 


0. 


013 





.059 


0. 


33 


0. 


73 


7.4 


UGC4879 


VV124 


dlrr 


MW 


1100 


7.3 




















Pegasus 


DD0216 


dlrr 


M31 


410 


7.2 


0. 


057 





.095 


0. 


10 


0. 


64 


7.4 


UGCA92 


EGB.0427+63 


dlrr 


MW 


1300 


7.2 




















Sag DIG 


ES0594-4 


dlrr 


MW 


1060 


7.1 


0. 


11 


0. 


.17 


0. 


20 


0. 


20 


11.6 


AndVII 


Cassiopia dSph 


dSph 


M31 


216 


7.1 


0. 


016 





.022 


0. 


022 


0. 


025 


12.9 


Andl 




dSph 


M31 


48 


7.1 


0. 


0038 





.0098 


0. 


087 


0. 


67 


8.9 


Andll 




dSph 


M31 


160 


7.0 


0. 


0049 


0. 


.0086 


0. 


076 


0. 


50 


9.2 


AndVI 


Pegasus dSph 


dSph 


M31 


266 


6.9 


0. 


0023 





.023 


0. 


19 


0. 


60 


9.0 


Leo A 


DD069 


dlrr 


MW 


800 


6.8 


0. 


13 





.31 


0. 


.65 


0. 


78 


6.2 


Antlia 




dSph 


MW 


1330 


6.8 


0. 


043 










0. 


43 


9.0 


Lcol 


DD074 


dSph 


MW 


270 


6.8 


0. 


0099 


0. 


.18 


0. 


50 


0. 


76 


6.4 


Aquarius 


DDO210 


dlrr 


MW 


950 


6.7 


0. 


.037 





.083 


0. 


12 


0. 


12 


12.0 


Andlll 




dSph 


M31 


68 


6.5 


0. 


0022 


0. 


.0061 


0. 


10 


0. 


47 


9.8 


Cctus 




dSph 


M31 


680 


6.4 


0. 


0045 


0. 


.013 


0. 


17 


0. 


52 


9.9 


LGS3 


Pisces 


dlrr 


M31 


284 


6.3 


0. 


.015 





.046 


0. 


16 


0. 


.43 


9.8 


LeoII 


DD093 


dSph 


MW 


205 


6.3 


0. 


0028 





.012 


0. 


025 


0. 


70 


8.8 


Phoenix 




dlrr 


MW 


405 


6.3 


0. 


027 


0. 


.071 


0. 


23 


0. 


42 


10.3 


Sculptor 




dSph 


MW 


88 


6.3 


0. 


.010 





.016 


0. 


026 


0. 


14 


12.6 


Tucana 




dSph 


MW 


870 


6.2 


0. 


0048 


0. 


Oil 


0. 


014 


0. 


30 


11.6 


AndXV 




dSph 


M31 


170 


6.2 




















AndXVI 




dSph 


M31 


270 


6.1 




















Sextans 




dSph 


MW 


86 


6.1 


0. 


00 





.00 


0. 


00 


0. 


00 


12.0 


AndV 




dSph 


M31 


117 


6.0 


0. 


0045 





.048 


0. 


066 


0. 


35 


10.8 


Carina 




dSph 


MW 


94 


6.0 


0. 


0065 


0. 


.0077 


0. 


43 


0. 


67 


7.1 


Draco 


DDO208 


dSph 


MW 


79 


5.9 


0. 


0004 


0. 


.010 


0. 


025 


0. 


19 


10.9 


Ursa Minor 


DDQ199 


dSph 


MW 


69 


5.9 


0. 


00 


0. 


.00 


0. 


00 


0. 


00 


12.0 


AndX 




dSph 


M31 


110 


5.9 




















AndXIV 




dSph 


M31 


162 


5.8 




















AndXVII 




dSph 


M31 


44 


5.8 




















AndIX 




dSph 


M31 


45 


5.7 





















Notes. — Listed are all presently known satellite galaxies of the MW and M31 within 1 h 1 Mpc of either host, with stellar mass 
M* > 5 x 10 J Mq . Star formation data are shown if available - here we parameterize the SFH in terms of the fraction of total 
stellar mass formed in the last 1, 2, 5, and 10 Gyr (i.e., /ig, /2G, /sg, /iog) and the mea n mass- weighted stella r age, r. Most of 
the star formation data are derived from the Local Group Stellar Photometry Ar c hive (Holtzman et al. 2006). Exceptions are 
M32, M33, Antlia, Sextans and Ursa Minor, which are taken fromlDolphin et al.l ((20051) and r is estimated f rom the reported 
/iog value. The LMC data are inferred from Smccker-Hanc et al. (2002); SMC from Harris & Zaritsky (2004). Distances from 
the host, rhost, are dire ctly taken or inferred from[Grcbel et al. (2003) with so me exceptions: A ndX ijZucker e t al. 200 7]), AndXIV 
llMaiewski et al .1 [20071') . AndXV and AndXVI dlbata et al.|[2007l) . AndXVII l|Irwin et al.ll2008ft, UGC487 9 (iKopvlov et al.ll20"08h . 
UGCA 92 (|Matedll998t ). and M33 (|McConnachie et alJ^Ool T Stellar mass, M„ is fr omlDekel fc Wod (|2003l ). except for M32 
and Sagittarius, whose values are estimated using M*/Lv = 3 Mq/Lq, with Ly from iMated (|l998); similarly we estimate M, 
for AndIX, AndX, AndXIV, AndXV, AndXVI, AndXVII and UGC 4879 from the quoted L v or M v in the same references 
cited for r hos t. 



4 



Orban et al. 




KGK04 (P KS = l.le-19) 
---e, = 0.1 (P K s = 0.0030) 
Local Group data 



1 



0.8 



^0.6 

V 

^0.4 



0.2 



I 1 


/" -J 








JT^ KGK04 (P KS = 2.0e-16)~ 




---e, =0.1 (Pats -0.43) 


1 


| Local Group data 



10"' 



/lG 



10^ 10" 

/2G 





KGK04 (P KS = 0.13) 
e. = 0.1 (P K s = 0.044) 
Local Group data 



0.4 0.6 

/lOG 



0.6 



0.4 



0.2 



KG^04 (P K s = O.OOlV) 




--e. =0.1 (Pk S = 0.61) 


/ ( ' r - ^ 


Local Group data 




- IV I 




i — I / > 




r f!~~ 




— f-f-~~---~r i i 





r (Gyr) 



10 



12 



14 



25 



15 





' ' KCKill O.'oOoij" 




---(,= 0.1 (F,rs = 0.0034) 




Local Group data 

y\ - 


.."I 


1 1 iii i^T"!-"" ~ 



10 



10 



10° 10° 

M, (Mo) 
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Fig. 1. — Top four panels: Cumulative distributions of the fractions of stellar mass formed in the last 1, 2, 5, and 10 Gyr. The f\a and 
f2G fractions reflect recent star formation, while the /5G and /iog fractions represent the overall star formation history. Solid line shows 
the data for the Local Group. Gray shaded region shows the spread of predictions of the KGK04 model, for 10 random realizations of the 
model. Dashed region shows the range predicted by our model with a stochastic star formation threshold, e» = 0.1, also for 10 realizations. 
The numbers in parentheses show the Kolmogorov-Smirnov probability of the model average being consistent with the data. Bottom left: 
cumulative distribution of the mass-weighted mean age of stellar population. Bottom right: Cumulative stellar mass function per host halo. 
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Fig. 2. — Three cumulative radial plots and a cumulative stellar mass function comparing the data for observed Local Group dwarfs 
(solid line) to three variants of the star formation model. The KGK04 model (gray shaded region) does not include stochasticity in the star 
formation density threshold, £ t hi (i- e -> e * = 0) whereas the two other models include stochasticity with a logarithmic dispersion e* =0.1 
(contours of dotted lines) and e* = 0.2 (contours of dashed lines). The other parameters of the Kennicutt-Schmidt law are kept the same 
for the three models. The models include an intrinsic spread from the random assignment of the angular momentum spin parameters to 
the halos. This spread is shown here as a filled region, or contours outlining a region, which encompasses the range of model predictions 
for 10 random realizations of the model. 



TABLE 2 

Stochastic Threshold of Star Formation Law 







e* = 


e, = 0.05 


e» = 0.1 


e, = 0.15 


e* = 0.2 


Pks- 


fall 


3.1c-l 


3.0e-l 


2.5e-l 


2.1e-l 


2.2c-l 


Pks- 


rdSph 


3.7e-l 


3.6e-l 


3.8e-l 


3.6e-l 


4.9e-l 


Pks- 


rdlrr 


8.7e-2 


5.2e-2 


4.2e-2 


3.3e-2 


3.9e-2 


Pks- 


flG 


1.1c- 19 


2.0c-6 


3.0e-3 


3.1c-2 


l.le-1 


Pks- 


Hg 


2.0e-16 


3.4e-2 


4.3e-l 


3.6e-l 


1.0e-l 


Pks- 


f5G 


3.8e-13 


2.5e-4 


2.1e-2 


7.8e-2 


2.9e-l 


Pks- 


/lOG 


1.3e-l 


9.9e-2 


4.4e-2 


1.6e-2 


2.5e-3 


Pks- 


T 


1.7e-3 


1.3c-l 


6.1e-l 


3.8e-l 


1.5c-l 


Pks-- 


Mt 


6.2e-3 


3.1e-3 


3.4e-3 


2.8e-3 


1.2e-3 


Dwarfs per halo 


20 


21 


22 


24 


26 



Notes. — Other fixed parameters are E t ho = 5 Mq pc 2 , / s f r = 1. The observed 
number in the Local Group is 46/2 = 23 dwarfs per host halo. 
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The last factor accounts for the less efficient gas dissi- 
pation in small halos with the virial temperature T vlr < 
few x 10 4 K, or equivalently, maximum circular veloc- 
ity V max < 50 km s _1 . This factor is written in terms of 
V4 = 16.7 km s _1 , the virial velocity corresponding to the 
virial temperature T v ; r = 10 4 K, and a constant, c = 10, 
chosen to reproduce the correct total number of dwarfs. 
This important factor, coupled with the density thresh- 
old for star formation, suppresses star formation in most 
low-mass halos. Equation @ is consistent with results 
of recent hydrodynamic simulations of galaxy formation 
(|Kravtsov k Gnedinl [20051 iTassis et all 120081 ) . An al- 
ternative parametrization of dwarf gas eous disks with a 
tempe rature floor at ~ 10 4 K is given bv lKaufmann et ahl 
(2007). In their model the disks are vertically puffed- up, 
which works to the same effect to reduce the gas density. 
The observed stellar core radii in the Local Group dwarfs 
are a factor of 2 — 3 smaller than those predicted by our 
model, which may favor the vertical expansion over the 
radial one assumed here. However, since it would not 
change the surface density S 9 , we cannot fit it in our 
framework of the Kennicutt-Schmidt law of star forma- 
tion. 

The spin parameter A is drawn randomly from a stan- 
dard log-normal distribution with A = 0.045 and <j\ = 
0.56 (jVitvitska et al.ll2002t [Hernandez et al.ll2007h . This 
adds an intrinsic variance in the predictions of the star 
formation model. While the KGK04 model used only 
one set of the randomly-selected A values to compare 
with the observations, in this paper we take 10 random 
realizations of each model, in order to account for this 
intrinsic variance. 

The gaseous disk is modeled spatially by 50 radial 
zones. At each simulation output epoch (about every 
10 s yr), the newly accreted gas is added to these radial 
zones, according to the exponential profile of eq. JT]). In 
each zone, the rate of star formation is determined by 
the Kennicutt-Schmidt law: 

±. = 2.5 x 10- 4 / sfr ( Sg ) M kpc- 2 yr"\ 
VI M pc 2 J 

(3) 

wherever the gas density, S ff , exceeds the threshold, Eth 
(|Kennicutti I1998D . Standard parameters are n « 1.4, 
/ s f r = 1, and E t h = 5 M Q pc~ 2 . Most of the variants 
of the model discussed in this paper employ this star 
formation law but with an important addition, described 
in S|4] below. 

In addition to this quiescent mode of star formation, 
the model also allows a starburst mode prompted by 
strong tidal interactions with other halos (see section 6.1 
in KGK04). In this mode there is no density threshold 
and even galaxies with E s < Eth can form stars if a 
strong enough interaction occurs. 

Finally, the e ffect of stellar evo l ution is taken into ac- 
count following iPrieto k Gnedinl (|2006[ ). They find that 
some 40% of the initial stellar mass is lost to stellar wind 
and supernovae after 5—10 Gyr, give n the assumpti ons of 
the initial stellar mass function from lKroupal (|2001h . stel- 
lar re mnant masses as a fu n ction of initial stellar mass 
from IChernoff k Weinberg (1993), and main sequence 
lifetimes from iHurlev et all (2000). In the models inves- 
tigated here we simply assume that 40% of the initial stel- 
lar mass is lost between the time when the stars formed 



and the present day. We do not "recycle" the liberated 
gas back into the ISM, which would make it available to 
form more stars. As a result our estimate of the stellar 
masses are the lower limit of the true masses - a simple 
decrease of the total stellar mass, M*, for each dwarf by 
a factor of 0.6. Accordingly, the /[i,2,5,io]G an d t val- 
ues are unchanged by stellar evolution in our models and 
since the HST photometry for the Milky Way satellites 
is typically good enough to measure the main sequence 
stars we simply compare the observed /[1,2,5,101c an( i r 
values (Table [T]) to the same quantities for the simulated 
dwarfs without any corrections. 

We mark as luminous (having corrected M* for stellar 
evolution as just described) those satellite halos with a 
predicted stellar mass > 5 x 10 5 M Q at z = 0. Such 
a cutoff agrees with the observational limit of all Local 
Group dwarfs (Table QJ known until a few years ago, 
when the ultrafaint dwarfs were discovered. We defer 
the discussion of the model predictions for these low- 
mass objects until $6] 

As in KGK04, we do a rough morphological classifica- 
tion of dwarf galaxies as dSph or dlrr based on the ratio 
of the stellar rotation velocity to the velocity dispersion: 
Vjot/v < 3 for dSph and v ro t/a > 3 for dlrr. In the 
model, the rotation velocity is calculated as the circular 
velocity at the outer-most stellar radius, as it would be 
measured in observation, while the velocity dispersion is 
estimated from the amount of external tidal heating in 
strong tidal interactions with other halos. Such classifi- 
cation does not take into account the recent star forma- 
tion activity or the remaining gas content, and therefore, 
is only a crude indication of the observationally defined 
morphological type. Comparison of these model predic- 
tions to the data is, in fact, completely complementary 
to the comparison of the star formation histories. 

In the process of revising our models we have discov - 
ered that the gas densities E s in iKravtsov et al.1 (|2004D 
were underestimated by a factor of 2 due to an error in 
the code. When we quote the results for the KGK04 
model here, we use the corrected values. 

3.1. Discrepancies of KGK04 Model with Star 
Formation Data 

Despite significant successes in explaining the number 
and spatial distribution of the Local Group dwarfs, the 
KGK04 model did not predict a sufficient amount of re- 
cent star formation. Figure Q] shows that 95% of the 
dwarfs in that model have not formed any stars in the 
last 1 Gyr, in serious disagreement with data. 

To quantify the level of this disagreement, we use the 
Kolmogorov-Smirnov (KS) test for cumulative distribu- 
tion functions of the following parameters: the distance 
to the host (r a n), the distances for the dSph and dlrr 
galaxies separately (rasph and rain-), the star formation 
fractions (f 1G , f 2G , f 5G , /iog), the mean age (r), and 
the stellar mass (M„). Table [5] shows that the recent 
star formation fractions in the KGK04 model (see column 
e* = 0) have very low KS probabilities, below 10~ 12 . On 
the other hand, the /iog fraction, which is a more global 
measure of the overall star formation, is consistent with 
the data at the 13% level. Thus the star formation law 
used in the model is not necessarily at fault, but the ap- 
parent lack of recent star formation is a direct result of 
the fixed star formation threshold, as we will see below. 
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The mean mass- weighted age is also inconsistent with 
the data, but at a less significant level (Pks ~ 10~ 3 ). 
The mean age is a global, integral measure of the SFH, 
which combines the low-probability recent star formation 
with a higher-probability early star formation. The stel- 
lar age is systematically overestimated in the model, by 
a few Gyr. 

The total stellar mass is also overestimated by a factor 
of several for most dwarfs. The cumulative mass function 
is inconsistent with the data at a level similar to the age 
distribution, Pks ~ 10~ 3 . 

Note also that KGK04 considered principally the satel- 
lite dwarfs located within the virial radius of the host 
galaxy at r < 200 h^ 1 kpc. We extend our analysis to 
all dwarfs in the Local Group, even the potentially iso- 
lated ones lying outside the virial radius of either host. 
Figure shows their cumulative radial distribution out 
to 1000 /i" 1 kpc. The KGK04 model does well for the 
distribution of all dwarfs and the distribution of Irr/dlrr 
types separately. The model slightly overestimates the 
number of dSph/dE types outside 200 h" 1 kpc, but in 
all cases the KS probability is ~ 10% or higher, fully 
consistent with the observations. 

Our phenomenological model contains several param- 
eters that allow for some freedom in the outcome. Two 
of the most important parameters in the model are the 
threshold density Sth and the disk size parameter, c (eq. 
[5]) , both of which arc difficult to know, a priori, from the- 
ory. Since these parameters can significantly change the 
predicted total number of dwarfs, the observed number 
of (M* > 5 x 10 5 Af Q ) dwarfs provides a fairly stringent 
constraint on these values. For example, with c = 10 
and Sth = 5 M Q pc -2 , the predicted average number of 
dwarfs per host halo within 1 Mpc is 20.5. If we take 
c = 5 (with S t h = 5 M Q pc -2 ) this number increases to 
40.7, and if we take c = 15 the number of dwarfs drops 
to 14.3. Analogously, reducing the threshold density to 
4 and 3 M Q pc -2 (with c = 10) increases the average 
number to 21.7 and 22.3, respectively. None of these 
fixed threshold models, however, adequately reproduces 
the observed star formation histories. Therefore, we look 
for additional physical ingredients for our model. 

In our extension of the KGK04 model, we attempt to 
retain the correct radial distribution of the dwarf galax- 
ies, while improving the predictions for their star forma- 
tion histories. 

4. STOCHASTIC STAR FORMATION THRESHOLD 

We consider a number of modifications to the KGK04 
model. The most promising of the modifications is the 
introduction of a stochastic threshold to the star forma- 
tion law. 

At each output epoch through the course of the simu- 
lation, the threshold density for star formation is drawn 
from a log-normal distribution with a mean value Etho 
and a small dispersion e* ~ 0.1: 



P(Eth) cEth = 



1 



2ttc, 



■ exp 



(logSth ~ logSthp) 
2e 2 



Hi- 



, (4) 

As a fiducial mean value we take S t ho = 5 M pc , but 
we also vary that parameter in some runs. 

What does this stochastic threshold mean? And what, 
physically, is the threshold of star formation? We ap- 



ply it to the azimuthally-averaged surface density of gas. 
In nearby star forming regions, which we can study di- 
rectly with HST and Spitzer, stars form in dense molec- 
ular clouds (e.g.. iMcKee k Ostrikerll2007t ). Locally the 
density is high, but the azimuthal average at a particular 
distance r from the center may be either high (if molec- 
ular clouds are common at r) or low (if they are rare). 
Thus, by invoking the threshold E t h, we arc effectively 
parameterizing the fraction of molecular gas available for 
star formation. When we take the threshold to vary, we 
are thus accounting for stochastic star formation in iso- 
lated HII regi ons, such as those found by GALEX in 
nearby spirals (jThilker et alJl2007f ). 

Nearby dwarf galaxies often show a very high HI 
gas fraction, i.e. the ratio of gas mass to baryon 
(gas+stars) mass, up to 90% (jFisher k Tullvi 119751 : 
iGeha et all 120061 ; iLee et~aT]|2006ft . The current star for- 
mation rates in those galaxies are low, and therefore 
most of the gas is inert. Hy drodynamic simulations of 
iRobertson fc Kravtsovi (|2008f l. which treat the formation 
of molecular hydrogen in detail, also suggest that dwarf 
disks may contain large reservoirs of diffuse atomic gas 
that is unable to condense in molecular clouds and par- 
ticipate in star formation. This inefficiency of forming 
molecular clouds at low densities effectively results in 
the threshold density for star formation. 

4.1. A Major Improvement 

Figure Q] shows that the stochasticity greatly improves 
the agreement of the recent star formation fractions 
(/ig>/2G) /sg) with the Local Group data. The pre- 
dicted distributions are much closer to the observed ones 
than in the KGK04 model, and are, in fact, statistically 
consistent with each other at ~ 10% level. 

We quantify this effect by gradually increasing the 
amount of stochasticity, e* . Table [2] shows the KS test 
results for ranging from 0.05 to 0.2, which approx- 
imately corresponds to 10% to 60% variation in the 
threshold density Sth- Even a small amount of stochas- 
ticity, e* ~ 0.05, leads to a dramatic improvement of the 
/1G1/2GJ and /5G statistics. The probabilities increase 
rapidly with e* to a good fraction of a percent or more. 
The model with = 0.1 is already statistically consis- 
tent with the data. 

With the other parameters of the star formation law 
(Stho and / s fr) being fixed, increasing leads to better 
recent star formation parameters but worse early star 
formation parameters. The probability of the /log dis- 
tribution decreases from 13% to under 5% for e* =0.1 
and even below a percent for = 0.2. 

The mean age is most consistent with the data (Pks — 
0.61) for e* = 0.1. The total number of dwarfs also in- 
creases systematically with e* while the stellar mass func- 
tion and the radial distributions are not significantly af- 
fected by the variation of e*. Additionally, the e* = 0.15 
and e* = 0.2 models begin to overpredict the total num- 
ber of luminous dwarfs per host halo, implying that at 
these values the overall star formation becomes too ef- 
ficient. This overabundance can be seen in Table [2] or 
graphically in Fig. [2] Therefore, we take the case with 
e* = 0.1 as our fiducial model striking the best balance 
for all star formation statistics. We consider other vari- 
ants of the stochastic model in S}5l 
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Fig. 3. — Left: The star formation histories of two simulated dwarfs in the KGK04 model, coarsened into 1 Gyr bins so as to better 
resemble the time resolution of the data. The x-axis shows the Age (i.e. of the isochrone) with the present epoch at Gyr, going back to 
the Big Bang at 14 Gyr. Right: The gas density profiles for the two simulated galaxies that were used to calculate their star formation 
rates, at five epochs. The fixed star formation threshold is shown by the dashed horizontal line. 
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Fig. 4. — Left: The star formation histories of the two dwarfs from Fig. [3] hi the model with a stochastic star formation threshold, 
e* = 0.1 (dashed histograms and quoted parameters). For comparison, shaded histograms show the observed SFHs of NGC 3109 and IC 
1613, normalized to the total stellar mass cited in Table [l] Right: The gas density profiles at five epochs. The range of variation of the 
star formation threshold is shown by dashed horizontal lines. This range allows star formation to continue at late times. 
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4.2. Reasons for Success 

Figures [3] and 0] illustrate why the stochasticity is so 
successful. They show star formation histories and cor- 
responding gas density profiles for two representative 
model dwarfs. Figure [3] is for the case of a fixed den- 
sity threshold, E t h = 5 M Q pc~ 2 . Without the stochas- 
ticity, and in the absence of radial gas flows, the gas 
above the threshold is steadily converted into stars at 
early times (10 and 12 Gyr ago) at substantial rates, 
0.1 — 0.3 Mq yr" 1 . At these epochs the galaxies evolve 
effectively in isolation, while still growing by gas-rich hi- 
erarchical mergers. At later time, when these galaxies be- 
come satellites of the larger host galaxy, they are tidally 
truncated and no new gas is accreted. After the high- 
density gas supply is exhausted, the rest of the gas hovers 
just under the threshold, unable to form new stars. In 
the first dwarf, shown in the top panels in Figure [3j star 
formation almost completely halts 8 Gyr ago. The other 
dwarf, shown in the bottom panels, experiences several 
distinct episodes of star formation, the last one finish- 
ing 3 Gyr ago. In both cases there is effectively no star 
formation in the last 2 Gyr, fic ~ /2G ~ 0. 

If instead the threshold density varies in time, at later 
epochs the gas may find itself above the threshold and al- 
low more recent star formation. Figure H] shows the star 
formation histories for the same dwarfs but now with a 
variable threshold and using e„ = 0.1. There is more 
star formation overall and more star formation at later 
times. As a result, the stellar mass is higher by 20% to 
30% and mean stellar age is lower by about 1 Gyr. More 
importantly, several percent of all stars are formed in 
the last 2 Gyr. This stochastic enhancement of the star 
formation rate allows the model to reproduce the star 
formation episodes at late times when gas-rich mergers 
(which increase the available gas supply) and strong in- 
teractions (which can prompt starbursts) are relatively 
less common. 

We also show on Fig. [4] the observed star formation 
histories of two dwarf galaxies, NGC 3109 and IC 1613. 
We do not expect the model predictions to correspond 
in detail to the observed SFH features, as our modeling 
is necessarily statistical and is aimed at explaining not 
a specific galaxy's SFH but only an ensemble of SFHs 
of many Local Group dwarfs. Still, qualitatively, the 
agreement between the model and the data is good: in 
the case of NGC 3109 both show most stars being formed 
at early times with a small fraction in the last 4 Gyr, 
while in the case of IC 1613, both SFHs continue until 
the present in several distinct, extended episodes. 

4.3. Remaining Discrepancies 

Despite the impressive improvements, there still re- 
main discrepancies of the fiducial model with the ob- 
served data. 

First, there is still not quite enough very recent star 
formation, a problem which is quantified by the KS-test 
result for / 1G in Table [2] (Pks ~ 3 x 10~ 3 ) and is ap- 
parent in Fig. [T] The models with a higher amount of 
stochasticity, e* = 0.15 and e* = 0.2, achieve better 
agreement with the observed fic distribution but they 
are disfavored for skewing all other star formation statis- 
tics. Thus, in the fiducial model (e* = 0.1) over 25% of 
the dwarfs have less than 10~ 3 of their stars formed in 
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Note. — Other parameters are e„ = 0.1, E t ho = 5 Mq pc , 

the last 1 Gyr, compared with about 10% of such dwarfs 
in the observed sample. Note that the current observa- 
tions are sensitive to SFHs with large fractions of star 
formation at late r times. The very deep HST ACS imag- 
ing of Leo A by I Cole et aTl (|2007t ) have confirmed that 
majority of all star formation has occurred in the last 
half of the age of the universe. 

Second, the total stellar masses of the luminous dwarfs 
in the models tend to be significantly above the observed 
stellar masses of the Local Group dwarfs. This result is 
apparent for any value of e* between and 0.2 (see bot- 
tom right panel of Fig. This excess stellar mass in our 
models can perhaps be traced to the fact that we do not 
include thermal and ionizing feedback from young mas - 
sive stars (e.g., iDekel fc Silkjfl98l iDekel fc Wod 12003). 
which can disrupt molecular clouds and drive galactic 
outflows, thus reducing the available gas supply for star 
formation. 

Third, the model does not predict enough very early 
star formation. This is a generic problem with any hi- 
erarchical model in which the presently-massive satel- 
lites form and accrete onto the host late and, at early 
times, therefore, do not contain significant amounts of 
gas above the threshold. Our model predicts that some 
20% to 30% of the dwarf galaxies will have formed greater 
than 85% of their stellar mass in the last 10 Gyr (i.e. 
/iog > 0.85). In other words, there exists in the model a 
number of dwarfs with anomalously young stellar popula- 
tions, whereas, by contrast, the Local Group data in Ta- 
ble [1] does not show any dwarfs having formed more than 
86% of its stellar mass in the last 10 Gyr (the Saggitarius 
dSph has the highest fraction) - at least 14% of all stars 
in the observed dwarfs formed in the first 4 Gyr after the 
Big Bang. Some of the simulated dwarfs with very young 
stellar populations include objects which form the bulk 
of their stars in one, punctuated star formation event 
caused by a close tidal interaction, prompting a star- 
burst which happened to occur in the last 10 Gyr. But 
not all of the galaxies in this category form stars through 
the starburst mode; other dwarfs with very young stellar 
populations form all of their stars through the continuous 
mode with no starbursts at all. 

5. OTHER VARIANTS OF STAR FORMATION MODEL 

With the discrepancies discussed in H4.3I in mind, we 
have explored several variants of our star formation 
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model in order to check if relaxing other model assump- 
tions can improve the predicted star formation histories 
to the point where there is broad agreement with the 
observations. 

5.1. Gas Accretion within the Virial Radius 

One of these alternate models is motivated by the lack 
of enough recent star formation. In the KGK04 model 
and in the stochastic models shown in Figs. [T] and ac- 
cretion of new gas onto the dwarf halos is shut off when- 
ever the dwarf comes within the virial radius of the host 
halo. This is based on the expectation that satellite ha- 
los near their host would be tidally truncated and unable 
to capture new gas even if they happen to increase their 
dark matter mass in mergers with other satellite halos. 
However, becoming a satellite does not immediately stop 
the halo's star formation activity which can continue as 
long as the gas density remains above the threshold or if 
a starburst event is triggered. 

In the first of these alternative models, we lift this as- 
sumption and allow the accretion of new gas even within 
the virial radius of the host, at d < i? V ir- As a result, the 
agreement between the model and the data for the fia 
distribution improves as indicated by the KS-test results 
shown in Table[3]- from the rejection of the null hypoth- 
esis of the same-distribution at 0.3% significance level 
to rejection only at the 9% level, essentially statistically 
consistent. 

Qualitatively, the improvement is most dramatic for 
fiG since allowing the additional accretion of gas is more 
important at late times, and because the virial radius of 
the host halo will be larger than at earlier epochs; also 
the dwarf galaxies at late times are likely to be closer in. 
The extra gas is converted into extra stars and so the 
fiG value increases and the typical mass of the dwarfs 
increases. Generally, the extra accretion shifts the stel- 
lar age distribution to be younger. However, another 
predicted feature of our models is a tail of the age distri- 
bution at r < 5 Gyr. The extra accretion exaggerates 
this tail since it increases recent star formation. In con- 
trast, the youngest dwarf galaxy in the Local Group is 
Leo A, with r = 6.2 Gyr (though see the Leo A SFH of 
ICole et al.ll2007l which has a significantly younger result 
for r). This is fundamentally the same problem as the 
existence of simulated dwarfs with /i g > 0.85, and the 
fwG distribution similarly becomes less consistent with 
the data in this model. 

5.2. Modifications of the Schmidt Law 

Another important question to be addressed is whether 
the observed SFHs, given our model assumptions, fa- 
vor steepening of the Schmidt-Kennicutt law (£* cx £"). 
There is some evidence for the exponents n to become 
grad ually larger at low gas densities in dwarf galaxies 
(e.g.. lBoissier et al.ll2003l: iHever et al.ll2004t see also dis- 
cussion in iRobertson fc Kravtsovll2008f) . Since most of 
our dwarfs are close to the threshold, we investigate this 
effect simply by changing the exponent at all densities, 
setting it to n — 1.4, 2, or 3. For each of these values, the 
normalization / s f r is allowed to float such that the overall 
star formation efficiency is lowered until the point where 
the median of the age distribution, T mo( j, of the simu- 
lated dwarfs matches that of the Local Group galaxies, 
9.25 Gyr. 
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Note. — Other parameters are e, = 0.1, Etho = 5 M© pc 

Table [4] shows the required values of / s f r and the re- 
sults of KS-tests. The rather high significance results 
from the KS tests applied to the r-distribution are not 
particularly surprising since the normalization of the star 
formation law has been fine-tuned for this comparison. 
Qualitatively this exercise keeps the overall star forma- 
tion rate the same, i.e. the stellar mass function is basi- 
cally unchanged. And similarly the radial distributions 
stay roughly the same - this is important since one of the 
primary successes of the KGK04 model was reproduc- 
ing the observed radial distribution of the Local Group 
dwarfs, something that our models continue to do. Even- 
tually the same amount of gas above the density thresh- 
old is converted into stars, so the rate at which star for- 
mation proceeds does not alter the total number or the 
radial distribution of the simulated dwarfs. 

By and large, all of the same discrepancies with the 
data listed in §4.31 are apparent in these models with 
different n - too large stellar masses, not enough recent 
star formation, and too much overall star formation in 
the last 10 Gyr. Arguably the {n = 3, / sfr = 0.012} 
model is a best fit to the data for fia and /sg, however 
the fiducial {n = 1.4, / s f r = 1} model is still the best for 
fiG and /iog- Further, comparisons of the stellar mass 
and the r-distribution are inconclusive as well. 

5.3. Revised Epoch of Reionization and Cutoff Mass 

Another interesting test is to vary the epoch of reion- 
ization. Appendix B of KGK04 present s an analytic fit to 
the numerical results of lGnedinT (|2000h that parametrize 
the reionization epoch in terms of the redshift when cos- 
mic HII regions begin to overlap, z Q , and the redshift 
when the reionization is completed, z r . Varying these 
two parameters allows us to test the sensitivity of model 
predictions to the details of the evolution of the ionizing 
background radiation. Note that our parametrization 
of reionization here applies only to the local neighbor- 
hood of the Galaxy and may differ from the global cosmic 
reionization. 

Reionization affects the amount of gas accreted onto 
dwarf halos. Increasing extragalactic UV flux during 
reionization photoionizes the gas inside and outside dark 
matter halos and prevents the halos with shallow poten- 
tial wel l s from capturing new gas heated to ~ 10 4 K. In 
iGnedinl (|2000f ) the effect of decreasing the gas fraction of 
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TABLE 5 
Reionization Scenarios 







Early 


Extended 


Fiducial 


Late 






9 


7 


7 


6 






10 


10 


8 


10 


Pks 




2.3e-l 


2.9e-l 


2.5&-1 


3.2e-l 


Pks 


fdSph 


5.7e-l 


5.6e-l 


3.8e-l 


5.3&-1 


Pks 


»*dlrr 


4.2e-2 


4.7e-2 


4.2e-2 


6.0e-2 


Pks 


flG 


1.6&-2 


7.7e-3 


3.0e-3 


1.3e-2 


Pks 


f2G 


7.7c-l 


6.4e-l 


4.3e-l 


6.6e-l 


Pks 


fsG 


1.2&-1 


8.2e-2 


2.1e-2 


6.3e-2 


Pks 


/lOG 


4.8e-3 


9.9e-3 


4.4e-2 


1.5e-2 


Pks 


T 


6.3&-1 


6.9e-l 


6.1e-l 


6.7e-l 


Pks 


Ah 


1.6e-3 


1.3e-3 


3.4e-3 


2.4e-3 


Dwarfs per halo 


20 


19 


22 


18 



Notes. — z r refers to the redshift when reionization is com- 
pleted, while z refers to the redshift where cosmic HII regi ons 
begin to overlap (see Appendix B in lKravtsov et a l. 2004 for 

details). Other parameters: e* = 0.1, Etho = 5 pc~ 2 , 

U = 1- 

halos was parametrized as 

/ gas = /b (1 + M c /Af)- 3 , (5) 

where /b is the universal baryon fraction, M is the halo 
mass, and M c is the cut-off mass parameter. In linear 
theory for baryon perturbations, this parameter can be 
related to the filtering mass as M c w 0.26 Mf. The filter- 
ing mass Mf is an integral of a function of the intergalac- 
tic gas temperature over cosmic history and corresponds 
to the mass of a halo that loses 50% of its baryons as a 
result of exter nal photoheatin g. The values of the filter- 
ing mass from|Gnedm| (pOOl can be calcu l ated f or any 
z a and z r using eq. (Bl) in lKravtsov etai] (l200l . 

We consider four scenarios listed in Table [5] early 
{z r = 9, z a = 10}, extended {z r = 7, z a = 10}, and 
late {z r = 6, z a = 10} reionizations, in addition to the 
standard scenario {z r = 7, z a = 8}, which was used in 
KGK04 and in our fiducial model. Generally, the results 
are quite similar for all scenarios, likely because lumi- 
nous satellites are hosted by the relatively massive halos, 
in which the virial temperatures are above 10 4 K and the 
external heating of the gas does not affect its distribution 
significantly. Early reionization is slightly preferred for 
the fie $2G, and /sg distributions, but statistically all 
scenarios are consistent with each other. The /iog dis- 
tribution is still reproduced best by the fiducial model. 

Recent hydrodynamic simulations indicate that the 
filtering mass may overestimate the mass of the ha- 
los that lose 50% of their baryons, especially at low 
redshift. T h e mos t sophisticated ART simulations of 
iTassis et al.l ((2008) , including the effects of radiative 
transfer, give distributions of / gas at several redshifts, 
z = 3.3,4,5,7,8,9 (their Fig. 2) for the standard model 
of reionization. We have fit these distributions with the 
form of equation ([5]) and obtained best-fitting values of 
M c . These values are similar to our old values at z > 8 
but deviate from them at lower redshift by as much 
as an order of magnitude. We supplement these high- 
redshift fits with th e z = results of SPH simulations by 
ICrain et al.l (|2007f ) and iHoeft etaU (|2007l ). which both 
give M c ~ 2 x 10 h^ 1 M Q , about a factor of 5 smaller 
than our old value. The combined sets at z = and 



z > 3 can be fit by the following expression, accurate to 
better than 50%: 

M c fa 1.8 x 10 7 + 3 x 10 9 (1 + z)- 3 h~ l M . (6) 

Ideally, we would like to derive the cutoff mass evolution 
from several simulations covering the whole range of red- 
shifts, but at the moment it is the best we can assemble 
from the literature. 

We have run our model with the new expression ((6j) 
for M c and found that the corresponding changes in the 
predicted SFHs are not straightforwardly better or worse. 
The fie distribution improves a little, f-xG and fee are 
effectively unchanged, while the /iog distribution is sig- 
nificantly worse off and the anomalously young dwarfs 
still persist. The age distribution and stellar mass func- 
tion are also more discrepant with the data than in the 
fiducial model. While the lack of improvement with the 
new cutoff mass prescription is unpleasant or, at any rate, 
suprising, until we have a more robust estimate of M c 
confirmed by several groups, we keep our fiducial model 
untouched. 

5.4. Monotonically Variable Star Formation Threshold 

In addition to exploring the stochastically variable den- 
sity threshold, we have also investigated a threshold E t h 
that varies with redshift monotonically. Though, at the 
moment, we do not know of a convincing physical mo- 
tivation for such a global systematic variation, we have 
investigated this possibility in exploring the full range of 
predictions of our model. Here we used a simple param- 
eterization of the S t h redshift dependence, 

£ th (z) = £ th0 (l + z) Q , (7) 

where £ t h(z = 0) = Stho = 3 M Q pc^ 2 and S t h (z — 
9) = 5 M pc- 2 , which results in a = log (5/3) ps 0.22. 

Interestingly, the results of this model are very simi- 
lar to our fiducial model with the stochastic threshold. 
The comparison is shown in Table [5] The new model 
shares the same problems as the fiducial model: the stel- 
lar masses are still too large and the anomalously young 
dwarf problem is still present (if not slightly worse) . We 
do not consider any other variants of the monotonic 
threshold but note that if a compelling theoretical or 
observational motivation for such a variation appears in 
the future, an evolving threshold may become a viable 
model. 



TABLE 6 

Stochastic vs. Monotonic Threshold 
of Star Formation Law 







e* = 
SthO = 5 


e* = 0.1 
SthO = 5 


<E* = 


Pks 


rail 


3.1c-l 


2.5e-l 


3.0e-l 


Pks 


^dSph 


3.7e-l 


3.8e-l 


3.4e-l 


Pks 


f'dlrr 


8.7e-2 


4.2e-2 


3.7e-2 


Pks 


flG 


Lie- 19 


3.0e-3 


3.7e-2 


Pks 


f-ZG 


2.0e-16 


4.3e-l 


7.1e-l 


Pks 


fsG 


3.8e-13 


2.1e-2 


1.6e-2 


Pks 


/lOG 


1.3e-l 


4.4e-2 


2.7e-2 


Pks 


r 


1.7e-3 


6.1e-l 


5.0e-l 


Pks 


M* 


6.2e-3 


3.4e-3 


2.0e-3 


Dwarfs per halo 


20 


22 


23 
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Fig. 5. — A comparison of three variants of the fiducial model which reject dwarfs that have not formed any stars by t re j e ct = 0.5 Gyr 
(filled gray shape), t rc j cc t = 1.0 Gyr (dashed line), or t re j ect = 2.0 Gyr (dotted line) after the Big Bang. Local Group data are shown 
as a solid black line. Note that in order to match the observed number of dwarfs, the star formation density threshold is lowered in these 
models, see Table[7] The only noticeable improvements, relative to the fiducial model, are in the f\oG distribution and in the corresponding 
lack of the low-r tail in the age distribution. 
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TABLE 7 
Rejecting Late Beginners 







^reject — 

0.5 Gyr 


^reject — 

1 Gyr 


^reject — 

2 Gyr 


No 
rejection 


SthO 




1.87 


2.6 


3.5 


5.0 


Pks 


rail 


5.3e-l 


4.0e-l 


3.4e-l 


2.5e-l 


Pks 


^dSph 


1.7o-l 


3.7e-l 


5.1e-l 


3.8e-l 


Pks 


»*dlrr 


2.2e-l 


9.6e-2 


7.3e-2 


4.2e-2 


Pks 


flG 


1.6e-3 


3.3e-3 


4.9e-3 


3.0e-3 


Pks 


f2G 


2.5e-l 


5.0e-l 


5.6e-l 


4.3e-l 


Pks 


fsG 


5.0e-3 


7.3e-3 


1.4e-2 


2.1e-2 


Pks 


/lOG 


3.3e-l 


3.5e-l 


3.2e-l 


4.4e-2 


Pks 


T 


1.3e-l 


3.1e-l 


4.0e-l 


6.1e-l 


Pks 


M, 


2.8e-4 


2.7e-4 


2.9e-4 


3.4e-3 


Dwarfs per halo 


22 


22 


23 


22 



Notes. — Here we reject dwarfs that have not formed any 
stars by i re ject = 0.5 Gyr, 1 Gyr, or 2 Gyr after the Big Bang 
(redshifts z ~ 10, z ~ 6, and z ~ 3.3, respectively). Other 
parameters: e* = 0.1, / 8 f r = 1. 



5.5. Rejecting Galaxies with Delayed Star Formation 

Finally, in directly addressing the problematic issue of 
the simulated dwarfs with too young stellar populations, 
we consider models that reject dwarfs with delayed star 
formation, i.e. the dwarfs that have not formed any stars 
within the first few Gyr after the Big Bang. 

The reasoning for such an ad hoc cut is motivated 
by the uncertainty in the detailed effect of the early 
UV background on the gas content of low-mass halos. 
The satellites without early star formation in our model 
must have acquired a significant gas reservoir only at 
late times. In the hierarchical paradigm, these satellites 
have been built by mergers of smaller objects, each of 
which carried an even smaller amount of gas. In order 
to become a galaxy in our model, at some time the to- 
tal gas density distribution in the satellite must reach 
above the density threshold. However, if for any reason 
we underestimated the effect of gas loss from small halos 
at early times, then at that time the combined amount 
of gas would be overestimated and star formation should 
not take place. Our treatment of the gas heating during 
reionization is very approximate, which makes it likely 
that we could either underestimate or overestimate the 
gas loss effect for particular dwarfs. In general, star for- 
mation at very high redshift could proceed even at lower 
densities than assumed in our models, such that all lumi- 
nous dwarfs form at least a fraction of their stars before 
reionization (|Ricotti fc Gnedinll2005h . 

In the current concordance cosmology (|Komatsu et al.1 
2008), the epoch of complete reionization can range from 
z^l0toz^6, which corresponds to a range of times 
from 0.5 to 1 Gyr after the Big Bang. Given this uncer- 
tainty, and uncertainties in the details of gas outflows, we 
consider three models that reject the dwarfs without star 
formation in the first 0.5 Gyr, 1 Gyr, or 2 Gyr after the 
Big Bang. The last model with i ro jcct = 2 Gyr is the least 
restrictive and closest to the fiducial model (no rejection), 
which can be formally written as i re joct > 14 Gyr. In 
each of these "reject late beginners" models, the star for- 
mation density threshold Stho is lowered in order to raise 
the number of dwarfs back to the number observed in the 
Local Group. The details of these models are listed in 



'reject = °' 5 'V r 

'reject = 1 G y 

'reject = 2 G y r 

no rejection 

o Local Group data 




Fig. 6. — Mass-weighted mean stellar age vs. stellar mass at 
z = 0. We show the fiducial model and three variants, which 
reject simulated galaxies with delayed star formation, by 0.5, 1, 
and 2 Gyr, respectively. The model galaxies are binned by stellar 
mass and only the bin averages are shown for clarity. Vertical bars 
indicate the standard deviation of the sample in each bin. Observed 
values for the Local Group (Table [TJ are plotted by circles and fall 
in the same range as the model. 

Table [3 

Figure O shows the star formation parameters in the 
three models. All of the distributions are similar to 
each other, but they all deviate from the fiducial model 
(plotted in corresponding panels of Fig. [TJ) in one im- 
portant aspect. The /ioG-distribution lacks objects with 
/iog = 1, i-e. with no star formation in the first 4 Gyr 
after the Big Bang. Unfortunately, the time resolution 
of the observed SFHs does not allow us to discriminate 
among the three variants of the cut. But all of them pro- 
vide the needed fix: the /ioG-distribution is fully consis- 
tent with the data (P K s > 30%). 

Table [7] shows that the other probabilities remain 
roughly the same or even decrease relative to the fidu- 
cial model. If we are to accept any of the "reject late 
beginners" variants, we would prefer the least restrictive 
Reject = 2 Gyr model. Ideally, of course, we prefer to 
develop a better understanding of early star formation 
that would make this cut unnecessary. 

Figure [S] shows the stellar age as a function of stellar 
mass in the fiducial model and its three variants. Reject- 
ing any of the "late beginners" increases the stellar age 
by less than 1 Gyr, which is significantly smaller than the 
dispersion of the sample at all masses. A general trend, 
also largely overcome by the dispersion, is for the mean 
age to increase with mass until M* ~ 10 9 M and then 
to decrease at larger masses. The Local Group data are 
consistent with the younger stellar ages, and extended 
SFHs, for more massive galaxies. At the lowest-mass 
end, however, the observations show a number of very 
old objects that are still not present in our models, even 
with the strictest age cut. It is apparent from this and 
previous plots that our prescription for star formation in 
the first few Gyr of cosmic time still needs improvement. 

It is possible that the starburst mode of star forma- 
tion is causing very young stellar ages in small galaxies. 
We have checked, however, that the anomalously young 
dwarfs do not all have significant starbursts at late times 
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TABLE 8 

Low Mass Satellite Galaxies of MW and M31 



100 



Galaxy 



Host r host (kpc) My M, (M Q ) Refi 



Leo T 

Canes Venatici I 

AndXI 

AndXIH 

AndXII 

Hercules 

BootesI 

Ursa Major I 

LeoIV 

Canes Venatici II 

SDSSJ100+5730 

SDSSJ1329+2841 

Ursa Major II 

Coma Bercncics 

BootesII 

WillmanI 



MW 
MW 
M31 
M31 
M31 
MW 
MW 
MW 
MW 
MW 
MW 
MW 
MW 
MW 
MW 
MW 



417 
218 
103 
95 
116 
138 
62 
106 
158 
151 
83 
76 
32 
44 
60 
38 



-8.0 
-7.9 
-7.3 
-6.9 
-6.4 
-6.0 
-5.8 
-5.6 
-5.1 
-4.8 
-4.2 
-3.9 
-3.8 
-3.7 
-3.1 
-2.5 



4.0 x 

3.7 x 

2.1 x 
1.5 x 

9.3 x 

6.4 x 

5.4 x 

4.5 x 

2.8 x 

2.1 x 

1.2 x 

9.4 x 

8.5 x 
7.8 x 

4.5 x 

2.6 x 



10 r ' 
10 s 
10 5 
10 5 
10 4 
10 4 
10 4 
10 4 
10 4 
10 4 
10 4 
10 3 
10 3 
10 3 
10 3 
10 3 



1 

2 
3 
3 
3 
4 
5 
6 
4 
4 
7 
7 
8 
4 
9 
10 



Refer e nces.— (11 Ide Jong et all 120081; (21 iZucker et all 
2006bl. [Martin et al l 1 20081: (31 iMartin et al l 120061: (41 
Belokurov ct al. 2007; (5) Bclokurov ct al. 2006, Sicgcl 2006; 
(6llSimon fc Gehall2007l; (71 iLiu et~alll2008l; (81 Zuckc r et al] 
2006al ; (91 I Walsh et al.ll2007l ; (10) Will man et al.ll2006l . ISiegel 
2006T 

Notes. — M* is estimated assuming M,/Ly = 3M /L Q , 
with Lv determined from Mv quoted in the references. 
Satellites of MW are discovered by the SDSS, satellites of 
M31 by the MegaCam survey. Satellite-to-host radii, ri lost , 
for AndXI, AndXII, & AndXIII are derived using their 
projected distances from M31, assuming dM3i = 785 kpc 
l|McConnachie et all 120051 ). 

and are hosted by dark matter halos with a wide range 
of masses. In fact, galaxies with the highest fraction of 
stellar mass built in starbursts are typically old, with r 
between 10 and 12 Gyr. Starbursts are not very impor- 
tant overall - 91% of the dwarfs have less than 10% of 
their stellar mass formed in the starburst mode. Also, we 
have varied the two parameters describing the starburst 
mode (the minimum required tidal force and the fraction 
of gas converted into stars) and found that it has little 
effect on the SFH distributions. 

It should also be mentioned that there does seem to 
be a trend, at z = 0, for the anomalously young dwarf 
halos to be more massive than the other, more typi- 
cal, dwarf halos in the model - observationally speak- 
ing, these dwarfs would have higher dynamical mass-to- 
light ratios. Interestingly, at the end of reionization these 
anomalous dwarfs are hosted by halos with a very wide 
range of masses, so that there is significant dynamical 
evolution from the epoch of reionization until the present 
day. In other words, the younger dwarfs really do have 
rather "tumultuous" lives. 

6. PROJECTIONS FOR LOW-MASS DWARFS 

Given the amazing rate of recent discoveries of the pop- 
ulation of ultrafaint dwarfs in the SDSS and MegaCam 
surveys, there stands a challenge to predict the yet-to-be- 
observed star formation properties of these objects. Our 
models should in principle be able to predict the mean 
age and stellar masses for the low-mass dwarfs, however 
current predictions are not satisfactory. 

Figure [7] shows our fiducial model (with the stochas- 
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Fig. 7. — The stellar mass function per host halo, extended down 
to 5 X 10 3 Mq. All dwarfs in our fiducial model at d < 1 h~ x Mpc 
are shown in gray. The Local Group data for M t > 5 X 10 5 Mq are 
shown with a thick solid line, while the thick dashed line shows the 
number of presently-known low mass (M„ < 5 X 10 5 Mq) dwarfs 
(all the SDSS and MegaCam dwarfs listed in Table [8). The thin 
dashed line is a likely extrapolation of the mass function to all sky, 
to account for the incompleteness of current surv eys. Thin solid 
line is an estimate of the luminosity function by IKoposov et al.l 
(2008), assuming M,/L v = 3. 

ticity parameter e* = 0.1) extended to masses as low as 
ill* = 5 x 10 3 Mq. This model reproduces most closely 
the observed SFHs of the higher-mass dwarfs, but at 
ill* < 5 x 10 5 Mq, it predicts only a modest increase in 
the number of galaxies. The thick dashed line in Fig. [7] 
shows the currently-known number of ultrafaint dwarfs, 
listed in Table [8] The model predictions actually agree 
very well with the observed number and with the gentle 
slope of the mass function. 

The problem, though, is that the SDSS data release 
5, where the bulk of the new dwarfs have been discov- 
ered, covers only 20% of the sky. Therefore, we might 
expect the full sky to contain 5 times as many yet- to-be- 
discovered ultrafaint dwarfs (| Simon & Gcha 2007). For 
the MegaC am survey, the inco mpleteness factor is even 
larger, ~ 9 (|Martin et al.ll2006h . Extrapolation to all sky, 
by multiplying the observed number of dwarfs by these 
correction factors, is shown by the thin dashed line in 
Fig. [3 Our model falls well below this corrected mass 
f unction. 

IKoposov et a l. (2008) calculate the expected luminos- 
ity function of the faint dwarfs more accurately, by es- 
timating the maximum accessible volume of the survey. 
They conclude that the luminosity function should rise 
as dN/dM v oc lO - 1 ^, which for a fixed mass-to-light 
ratio results in N(> Af„) oc M ~ ' 25 . This estimate is also 
plotted by a straight line in Fig. [7] It lies below our first 
naive estimate but still significantly above the range of 
the model. Both incompleteness corrections predict over 
50 satellites above 10 4 Mq, a factor of two larger than in 
the model. Plus, more dwarfs may remain undetected at 
larger distances than those probed by the current surveys 
(~ 300 kpc). 

The discrepancy at low mass persists for all variants 
of our model and we were unable to find a set of pa- 
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rameters which adequately reproduced the stellar mass 
function at M* > 5 x 10 5 M© while predicting apprecia- 
ble numbers of M* < 5 x 10 5 M Q objects. This rather 
striking result is evidence that our models are not cap- 
turing some critical aspects of the formation of very low 
mass galaxies. The most naive solutions of significantly 
reducing either the threshold S t h or the disk structure 
factor c in eq. ([2]), thereby allowing a larger portion of 
gas to participate in star formation, are not viable op- 
tions since these modifications invariably overpredict the 
number of Af* > 5 x 10 5 M© dwarfs. In other words, 
lowering the density threshold produces so many dwarfs 
that the "missing satellites" are no longer missing and we 
are once again left with the expectation that we should 
see in the sky ~ 75 or more luminous dwarf galaxies 
around the Milky Way, whereas we only observe ~ 30. 
iKangj (|2008[ ) also reports a deficit of ultra-faint dwarfs 
from an independently-developed semi-analytical model 
of galaxy formation applied to a different collisionless TV- 
body simulation and using the same cosmic reionization 
model employed here. 

Observationally, a new interesting puzzle for our un- 
derstanding of dw a rf gal axy formation is presented by 
iRvan- Weber et all (|2008f) . In that study they observe 
the HI emission from Leo T, the only ultrafaint dwarf 
with measurable gas content and recent star formation. 
Based on the density and temperature of HI gas and ve- 
locity dispersion of stars, they find that the gas is every- 
where globally Jeans-stable, whereas the observed pock- 
ets of blue, 200 Myr-old stars indicate continuous star 
formation. The observed peak of HI column density is 
a few M©pc -2 , close to the star formation threshold. 
Leo T may thus present another example of stochastic 
star formation in a handful of isolated molecular clouds, 
surrounded by largely inert atomic gas. Another rela- 
tively massive dwarf, Canes Venatici I, also shows a small 
fraction of relatively young (~ 2 Gyr) , more metal-rich 
stars in addition to the predominantly old (~ 12 Gyr), 
metal-poor population (Ma rtin et al.l (2Q08h . Thus even 
the ultrafaint dwarfs may, in the future, reveal complex, 
extended star formation histories. 

7. CONCLUSIONS AND DISCUSSION 

We have presented phenomenological models for star 
formation histories of dwarf galaxies in the Local Group, 
based on the mass assembly histories in cosmologi- 
cal simulations and the stochastic density threshold in 
Kennicutt-Schmidt law of star formation. Our main con- 
clusions are as follows: 

• Models with a stochastic star formation threshold are 
much more successful than non-stochastic models, such 
as KGK04, in reproducing the observed star formation 
histories of the Local Group dwarfs. While the KGK04 
model predicted 95% of luminous dwarfs without any re- 
cent (t < 1 Gyr ago) star formation, our fiducial stochas- 
tic model correctly predicts that most dwarfs form a few 
percent of their stellar mass at late times, in agreement 
with the recent star formation fraction inferred for the 
Local Group dwarfs (see Table[IJ f\a column and Fig.[T]). 
Stochasticity allows star formation to proceed in isolated 
regions at late times if the threshold decreases. 

• Despite significant improvements of the stochas- 
tic model, a some discrepancies with the data remain. 
(1) Total stellar masses are typically too large by a fac- 



tor of several. This is a generic problem of both the 
fixed-threshold and stochastic models presented here. 
(2) About 10% of the dwarfs in both fixed-threshold and 
stochastic models have anomalously young stellar popu- 
lations. These objects build the bulk of their stellar mass 
in the last 10 Gyr, whereas all of the observed dwarfs 
contain at least 15% of stars older than 10 Gyr. These 
young stellar populations in the models are not created 
only by tidally-induced starbursts, but rather represent 
late mass assembly of some of the larger satellites. 

• Relaxing several model assumptions does not signif- 
icantly alter these predictions. We have considered the 
following variants of the fiducial model: allowing late gas 
accretion within the virial radius of the host halo; differ- 
ent slopes of the star formation law; extended epoch of 
reionization; and different prescriptions for the photoc- 
vaporation of gas from low-mass halos after reionization. 
All of these variants predict statistically similar observ- 
ables to the fiducial model. 

• A variant of the fiducial model that rejects dwarfs 
with no star formation in the first 1 or 2 Gyr after the 
Big Bang significantly improves the /iog stellar frac- 
tion. Even though we do not yet have an adequate jus- 
tification for such a cut, this model predicts the /iog- 
distribution fully consistent with the Local Group data. 
Stellar masses, however, are still overestimated. 

• Our fiducial model predicts only a modest popula- 
tion of dwarfs with M* < 10 5 M©, such as those re- 
cently discovered by SDSS and the MegaCam survey. 
The predicted stellar mass function would be an under- 
estimate if the observed numbers are extrapolated to all 
sky. However, our mass function is still consistent with 
the presently known dwarfs. 

Our phenomenological model contains several free pa- 
rameters, which allows significant freedom in the range of 
predicted properties of the satellite galaxies. As we dis- 
cuss in £)3.1l a combination of the parameters, c and S t ho, 
is well constrained by the observed number of dwarfs. 
Other parameters of the model, apart from the stochas- 
ticity e* which strongly affects late-time star formation, 
lead only to small and relatively insignificant varations 
from the fiducial model. 

This analytical prescription for star formation is ap- 
plied to the mass assembly histories of the halos in the 
cosmological A-body simulation. Since we average over 
three host halos, our results should not depend signifi- 
cantly on a particular host-halo merger history. Also, by 
construction all the massive satellite halos that become 
galaxies in our model survive tidal disruption in the host 
halo, so that their total number is predicted robustly 
even if their mass after tidal stripping may depend on 
the particulars of the simulation. 

We should emphasize that our model is not a unique in- 
terpretation of the SFH data. Our inferences here necce- 
sarily depend on the assumptions of a Schmidt law of star 
formation and an exponential profile for the gas density 
distribution inside dwarf halos as well as the assumptions 
we have made with the minimum density threshold for 
star formation, E t h- With this latter part of the model 
we have found that both stochasticity in this threshold 
over time (<2J) and a monotonically decreasing function 
f £|5.4[) can lead to extended SFHs - the feature missing 
from the KGK04 model. 
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We have assumed in our model that the gas clouds 
moving on circular orbits should generally remain at the 
same distance from the galaxy center. However, during 
mergers and tidal interactions the angular momentum of 
the gas can be perturbed, leading to radial infall. Such 
inflow of gas may bring the central gas density above 
the threshold in some halos soon after the Big Bang and 
lead to more early star formation. Modeling this process 
would be very interesting, as it may provide a nice so- 
lution for the /iog problem, but the complexity of such 
a process is beyond our simple model and requires a de- 
tailed hydrodynamic simulation. 

Additionally, we have not considered any gas outflows 
due to the radiative and thermal feedback of young stars. 
Such processes could reduce the amount of available gas 
supply and the total stellar mass of the simulated dwarfs, 
which may lead to a closer agreement of the predicted 
and observed stellar mass function. However, given the 
uncertainty in the evolution of the gas density profile, we 
cannot conclude that such feedback is required to recon- 
cile the stellar masses. In fact, from Figure [7] we see that 
in very low mass halos, where feedback is expected to 
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be stronger, our model needs a boost, rather than a re- 
duction, of star formation. The actual complex details 
of the condensation of molecular clouds in dwarf halos, 
which lead to the formation of stars, may turn out more 
important than the feedback of the stars after their for- 
mation. 

In spite of these uncertainties, the conclusion we would 
like to draw from our investigation is that the complex 
and extended SFHs of the Local Group dwarfs are gener- 
ally consistent with the expected star formation in cold 
dark matter halos, and that this star formation is gen- 
erally governed by the (low) efficiency of conversion of 
atomic gas into molecular clouds. 
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